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ABSTRACT 

We use 62,185 quasars from the Sloan Digital Sky Survey DR5 sample to explore the 
quasar mass-luminosity plane view of virial mass estimation. Previous work shows de- 
viations of ~ 0.4 dex between virial and reverberation masses. The decline in quasar 
number density for the highest-Eddington ratio quasars at each redshift provides an 
upper bound of between 0.13 and 0.29 dex for virial mass estimates. Across different 
redshift bins, the maximum possible Mgll mass uncertainties average 0.15 dex, while 
H/3 uncertainties average 0.21 dex and ClV uncertainties average 0.27 dex. Any phys- 
ical spread near the high-Eddington-ratio boundary will produce a more restrictive 
bound. A comparison of the sub-Eddington boundary slope using H/3 and Mgll masses 
finds better agreement with uncorrected Mgll masses than with recently proposed cor- 
rections. The best agreement for these bright objects is produced by a multiplicative 
correction by a factor of 1.19, smaller than the factor of 1.8 previously reported as 
producing the best agreement for the entire SDSS sample. 
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1 INTRODUCTION 

Strong correlations between the mass of supermassive 
black holes (SMBH) and the stellar velocity dispersion 
. l|Ferrarese fc Merritt 20001 IGebhardt et al. 2000|) and lumi- 
nosity (Ferraresc & Ford 2005J of their host galaxies ar- 
gue that our understanding of galactic formation is in- 
complete without an understanding of the SMBH found 
at their centres. Our modern arsenal for learning about 
SMBH growth is predicated on four basic tools: (1) Large 
samples containing ~ 10 5 quasars provided by mod- 
ern surveys (|Schneider et al. 20071 ISkrutskie et al. 20061 

ICroom etT al. 2004); (2) Bolometric luminosity estimation 
comparing a piece of the spectrum (Richards et al. 2006) 
to templates made from composite quasar spectra 
l|Elvis et al. 1994p : (3) The Soltan argument ||Soltan 19821 

ISalucci et al. 19991 lYu fc Tremaine 2002 [) that the inte- 
grated luminosity density of active galactic nuclei is consis- 
tent with the mass density in the local SMBH population; 
and (4) A series of SMBH mass estimation techniques com- 
prising the "BH Mass Ladder" (|Peterson et al. 2004)) . 

Until recently, the first three of these tools might 
have been thought sufficient. The highest-redshift rung of 
the black hole "ladder", virial masses, is the only op- 
tion at redshift z > 0.2, and comparison with rever- 
beration mapping yields a ~ 0.4 dex statistical uncer- 
tainty (Vestcrgaard & Peterson 2006). Further, early con- 



structions of the quasar mass-luminosity distribution showed 
a population of quasars all within ~ 1 dex of their Eddington 
luminosity ( Kollm eier et al. 2 006 ) . 

If quasars are all within ~ 1 dex of their Eddington 
luminosity (Kollmeier et al. 2006), then quasar luminosity 
would be a good proxy for supermassive black hole (SMBH) 
mass. Quasars could then be modelled as "light-bulbs" of a 
characteristic wattage, either operating near their Edding- 
ton luminosity or lying dormant (although Hopkins & Hern- 
quist 2008 argue this is an oversimplification). If the rela- 
tionship between quasar mass and luminosity were truly this 
simple, virial mass estimation would have little to add to the 
information already contained in quasar luminosity func- 
tions predicated on higher-precision bolometric luminosity 
estimates. 

In Papers I and II (|Steinhardt fc Elvis 20091 
ISteinhardt fc Elvis 2 009 ). we discovered that the quasar 
mass-luminosity plane is quite complex. We demonstrated 
the existence of several surprising new boundaries in the 
mass-luminosity plane. So, quasar mass and luminosity 
functions are not enough; these projections of the quasar 
mass-luminosity plane hide information. In Paper II , we 
showed that the ratio of emission lines from the broad 
line region may be evolving as quasars approach turnoff. 
High-precision black hole mass estimation would allow 
the identification of specific quasars very close to these 
boundaries as part of a search for a signature of their 
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underlying physical causes. However, SMBH masses cannot 
simply be well-estimated from the quasar luminosity, and 
so the investigation of SMBH evolution indeed requires all 
four tools. 

Our initial work relied upon H/3- and ClV-based virial 
mass estimates from Vestergaard & Peterson (2006) and 
Mgll-based estimators from McLure & Dunlop (2004). 
There have been several attempts to improve these 
virial mass scaling relationships (Onkcn & Kollmcicr 2008; 
|Risaliti, Young, fc Elvis 2009| IMarconi et al. 2009|l . In this 
paper, we reconsider virial mass scaling relations in light of 
these new boundaries in the quasar M — L plane. 

Every boundary in the SDSS quasar locus in the mass- 
luminosity plane is a combination of a real, underlying phys- 
ical boundary, SDSS selection, statistical uncertainty, and 
systematic uncertainty. The Eddington luminosity and sub- 
Eddington boundary (SEB) characterise the quasars with 
the highest Eddington ratios at each mass, or equivalently, 
quasars with the highest Eddington ratio at fixed luminos- 
ity. While we have termed this boundary "sub-Eddington" 
because of its behaviour at high mass, the SEB is flatter 
than the Eddington limit and for the lowest-mass objects 
at most redshift, the Eddington luminosity is the more re- 
strictive bound. As discussed in Paper I, the SDSS selection 
function is not a factor in the location of the SEB. Therefore, 
the decline in number density is a combination of statisti- 
cal uncertainty, systematic uncertainty, and an underlying 
physical cutoff, and as such, its sharpness can be used to 
derive an upper bound on the maximum statistical uncer- 
tainty in virial mass estimation. The statistical uncertainty 
has been estimated as 0.4 dex via comparison of virial and 
reverberation masses (Vestergaard & Peterson 2006). 

In § [2] we use the sharpness of boundaries in mass to 
place a tighter upper bound on the statistical uncertainty of 
virial mass estimates. In §[3l we show that proposed correc- 
tions to Mgll-based virial mass estimation are inconsistent 
for the brightest quasars at 0.6 < z < 0.8. We discuss the 
implications of these results in § [4] 



2 A NEW BOUND ON STATISTICAL 
UNCERTAINTIES IN SMBH MASSES 

We use black hole masses for 62,185 of the 77,429 SDSS DR5 
quasars (Schneider ct al. 2007) as determined by Shen et al. 
(2008) using H/3- and ClV-based virial mass estimators from 
Vestergaard & Peterson (2006) and Mgll-based estimators 
from McLure & Dunlop (2004). We divide the SDSS quasar 
population into 14 redshift bins, of width 0.2 below z = 2 
and wider at higher redshift. In each bin, we consider the 
mass distribution of quasars in a luminosity bin of width 
0.2 dex centred at the average bolometric luminosity for the 
catalogue at that redshift. Choosing this fixed luminosity 
gives us a large number of objects in our attempt to fit the 
decay rate at low mass. The decay at low mass is more rapid 
than at high mass. In each bin, we fit the decline in number 
density as both a Gaussian and an exponential decay (Figure 
[1} , reporting the dispersion or e-folding in Table Q] 

What form should we expect from the decay in quasar 
number density? Let the true quasar mass distribution be 
p p (M) and let 4>(x) be the probability distribution that a 
virial mass is incorrect by x. Then, the observed quasar mass 
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Figure 1. Top: The quasar M — L plane for Mgll masses at 
0.4 < z < 0.6. We take the quasar number density as a function 
of mass around the average luminosity (between the red, dashed 
lines). Quasar luminosity variability of 0.3 dex might lead to a 
change in estimated virial mass of 0.15 dex (blue). Bottom: The 
low-mass decay in quasar number density as a function of mass at 
45.46 < logLf, oi < 45.66 is best-fitting by an exponential decay 
with an e-folding of 0.15 dex (red). For comparison, the best- 
fitting decays using the average e-foldings for H/3 (cyan, 0.21 dex) 
and CIV (magenta, 0.27 dex) are also indicated. 



Table 1. Best-fitting forms for the decay in quasar number den- 
sity as a function of mass. These e-foldings imply maximum sta- 
tistical uncertainties of, on average, 0.21 dex for H/3-based virial 
masses, 0.15 dex for Mgll, and 0.27 dex for CIV. 
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distribution p {M) is the convolution 

/•CO 

p (M) = ( Pp * 0)(Af) = / p P (p)4>(M - p)dp. (1) 



For example, if p p (M) were a step (Heaviside) function, i.e., 
a constant number density at low mass and zero above the 
SEB, and the uncertainty <j>{x) were Gaussian, we would see 
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a decay taking the form of the error function, 1 + Erf(x) — 
e~* dt. In practice, the exact form of the tail is highly 
sensitive to p p , and the exponential decay described in Table 
[T]is a better fit to the low-mass tail of p than a Gaussian 
(Table [TJ, polynomial, or Erf(x). The exact form of p p ap- 
pears to be more complicated than a step function at masses 
above the SEB (Figure [TJ. However, the convolution acts to 
spread out the signal, and therefore the dispersion of fea- 
tures in p will be at least as large as those in <j>(x). So, the 
e-foldings of best-fitting exponential decays in Table [JJ are 
inconsistent with an 0.4 dex statistical uncertainty for virial 
mass estimation. 



2.1 Effects of quasar variability on virial mass 
estimation 



Virial mass estimates take the form 



log(M/M Q ) = A 
+ log 



/ FWHM(BLR line) \ 2 / XL X (B A) 
V 1000 km/s J V !0 44 erg/s 



for emission lines in the broad line region (BLR) and a 
nearby continuum flux. For the virial mass estimates used 
in the Shen et al. (2008) catalogue, C is between 0.50 
and 0.53 QMcLure & Jarvis 20021 |McLure fc Dunlop 2004] 
|Vestergaard fc Peterson 2006[ ). One possible explanation for 
this discrepancy, then, is that the analysis above only con- 
siders quasars at fixed bolometric luminosity, while the adja- 
cent continuum is one component of virial mass estimation. 
Perhaps the statistical uncertainty in virial mass estimation 
is mostly caused by quasar variability, in which the bolo- 
metric luminosity (and adjacent continuum) change on a 
timescale of years while the black hole mass remains very 
nearly constant. 

Variability moves objects along a L = M 2 line (Figure 
[TJ, blurring each thin luminosity slice along the mass axis 
by 0.18 dex for typical long-term, 0.3 dex optical variations 
(|de Vries et al. 2005[) . This blurring is larger than most of 
the Mgll decays, notably smaller than the CIV decays, but 
is consistent with the H/3 decays. The additional CIV scat- 
ter must have some other cause. Variability is wavelength 
dependent in quasars, being stronger toward the UV, but 
the difference between 2800 A and 1500 A is too small to 
explain the different Mgll and CIV decays 

Similarly, the adjacent luminosity used to estimate the 
black hole mass is based on the continuum local to each line, 
while the bolometric luminosity is based on the five SDSS 
photometry points. Thus, systematic offsets between these 
two luminosities could cause a scatter. The correction would 
have to have a peculiar shape as the longest and shortest 
H/3 and CIV measurements both have larger scatter than 
the intermediate wavelength Mgll value. 



3 THE MGII VIEW OF THE 

SUB-EDDINGTON BOUNDARY 

Mgll has the sharpest boundary in the M — L plane, and 
so appears to be the most precise mass indicator. Currently, 
though, H/3-based masses are considered to be the most re- 
liable, because they have been calibrated directly against 



reverberation masses at low redshift. Certainly, several po- 
tential problems with CIV virial masses have been suggested 
(|Shen et al. 20081 [Marconi et al. 2009)) . 

Corrections have also been suggested to the Mgll masses 
derived by McLure & Dunlop (2004) (MD04). Onken & 
Kollmeier (2008) examine SMBH for which both H/3 and 
Mgll masses are available and find that the Mgll-based Mbh 
may be overestimated at high Eddington ratio and under- 
estimated at low Eddington ratio. Risaliti, Young, & Elvis 
(2009) (RYE09) quantify this correction empirically as 



log[M BH (Hf3)} = 1.8 x ]og[M BH (MgII)] - 6.8. (3) 



In addition to correcting the central values of Mgll virial 
(2)masses, the RYE09 correction also increases the e-folding 
decay by the same factor of 1.8, to an average of 0.27 dex, 
which would make Mgll masses less precise than H/3 and 
comparable to CIV. 

This multiplicative correction is surprising, because it 
requires that the gas emitting either H/3 or Mgll lines 
has a non-virial component. If H/3 is virial, as expected 
from calibration between H/3 and reverberation masses 
( |Vestergaard fc Pete rson 2006), then the mass- velocity re- 
lation for Mgll would need to be M oc v 3 6 . RYE09 propose 
that this mismatch might instead be due to uncertainties in 
measurement of the Mgll line, mainly because of potential 
Fell contamination. 

To test the multiplicative correction we examine the 
SEB more closely. The slope of the SEB is sensitive to the 
multiplicative correction to Mgll-based virial masses but not 
to any additive correction. As the SEB is composed of the 
brightest quasars in each redshift bin they typically have 
the best-measured spectra at each mass, minimising Mgll 
measurement errors. Both Mgll and H/3 masses are available 
at the 0.4 < z < 0.8. Using the techniques detailed in Paper 
I we subdivided this bin into four redshift bins of width 
Az= 0.1. We then calculated the slope of the best-fitting 
SEB in each bin using both the MD04 scaling relation and 
the RYE09 correction. In Figure [2] we show the four H/3 
SEB slopes (green) compared to slopes using MD04 (blue) 
and RYE09 (magenta). 

In each redshift bin, the deviations between MD04 and 
H/3 slopes are between 0.3 a and 1.6 a (MD04 averages 0.8cr 
higher). Deviations between RYE09 and H/3 are between 
0.9 and 2.7 a (average 1.9 a), with the RYE09 slope always 
lower. If Mgll masses with the RYE08 correction produced 
identical SEB slope estimates to H/3, the probability that all 
four measurements would be this far below the H/3 slope is 
0.2% (comparing to a Monte Carlo of normally-distributed 
measurements). A best-fitting correction between H/3 and 
Mgll slopes (Figure [2]) reduces the factor of 1.8 to one of 
1.19. The M — L plane view of Mgll virial masses for the 
brightest quasars at each redshift is that corrections might 
not be necessary, and if necessary, are likely substantially 
smaller than previously proposed. Most of the RYE09 cor- 
rection may indeed be a result of measurement difficulties for 
Mgll, which become unimportant for the brightest quasars 
with the best-measured spectra. 
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Figure 2. A comparison of best-fitting slopes to the sub- 
Eddington boundary in four redshift ranges using different virial 
mass scaling relations. The H/3-based estimates (green) are a bet- 
ter fit for the original McLure & Dunlop (2004) Mgll masses 
(blue) than the Risaliti, Young, & Elvis (2009) correction (ma- 
genta). The best match between slopes is produced by a smaller 
Mgll correction (red). The black dashed line is drawn at the slope 
of the Eddington luminosity 



4 DISCUSSION 

The rapid falloff in quasar number density at fixed lumi- 
nosity near the sub-Eddington boundary (SEB), places a 
strong upper bound on the statistical uncertainty of virial 
mass estimation, of, on average, 0.21 dex for H/3-based virial 
masses, 0.15 dex for Mgll, and 0.27 dex for CIV. We consid- 
ered the possibility that the narrow spreads were induced by 
quasar variability. We find that variability can explain only 
~0.15 dex of spread. 

The small scatter in Mgll-based virial masses at the 
SEB implies, surprisingly, that the Mgll masses are more 
reliable than the H/3-based masses. We also find that, at least 
for the most luminous objects in each redshift bin, any Mgll 
to H/3 mass corrections are likely smaller than previously 
thought. 

We note that the ordering of the decay sizes is in- 
verse to the ordering of the emission lines by distance from 
the ionising continuum, as derived from reverberation map- 
ping (Peterson ct al. 2004). This ordering suggests a phys- 
ical explanation in which the non-virial motions in the 
broad-line region become relatively weaker with increas- 
ing distance from the continuum. An obvious candidate 
is a reduced role of radiation pressure at larger distances 
(|Marconi et al. 2009P , which would affect CIV most strongly 
and Mgll the least. 

The measured decays provide a fixed error budget to 
be distributed between non-virial motion and quasar lumi- 
nosity variability. If fluctuations in the adjacent continuum 
scale linearly with the quasar bolometric luminosity, the 
~ 0.15 dex bound on Mgll masses limits quasar variability 
to ~ 0.33 dex. If the decay rate is due primarily to non-virial 



motion in the broad line region, then the contribution from 
quasar luminosity variability is smaller. 

However, if the statistical uncertainty is truly smaller in 
virial masses than previously believed, then we must explain 
the origin of the ~ 0.4 dex scatter between virial and rever- 
beration masses, which would have to be dominated by un- 
certainties in reverberation masses. This situation is not easy 
to understand. It would be remarkable if the continuum fit at 
one wavelength from a single-epoch spectrum (virial masses) 
provides a more precise radius indicator than a time-series of 
high-resolution spectra (reverberation masses). Both virial 
mass estimation and reverberation mapping are based on as- 
sumptions about quasar geometry, so some of these assump- 
tions would need to be wrong. The basic difference between 
reverberation and virial masses lies in how each method de- 
termines the radius to broad emission lines, so the larger 
statistical uncertainty may lie in using time delays to esti- 
mate the radius. 

Reverberation data is also used to calibrate the rela- 
tion that allows the use of adjacent continuum as a proxy 
for radius in virial masses. A restriction to the best reverber- 
ation data includes scatter corresponding to typical uncer- 
tainties of 0.09 dex for broad line region radii from reverber- 
ation mapping, with an 0.11 dex uncertainty for the H/3 line 
( |Peterson, B. M., private communication) . If the scatter in 
the R — L relation is uncorrelated with other errors, this 
would restrict the remaining uncertainty to just ~ 0.1 dex 
for Mgll. Further, these uncertainty measurements use an 
improved R — L relation ( Bcntz et al 2009 ) , while the Shen 
et al. (2008) catalogue does not include these improvements. 

The two methods also differ in that virial estimates use 
the entire emission line, while the best reverberation meth- 
ods use the rms line profile, i.e., the minority fraction of the 
emission line that responds to continuum changes. Korista et 
al. (2001) show that only the gas that is optimally emitting 
a given line will respond to continuum changes. The virial 
method averages over gas in the entire broad line region. 
In some manner, virial masses appear to be giving better 
values than reverberation masses. We cannot satisfactorily 
explain this surprising result. 



5 CONCLUSIONS 

We used the sub-Eddington boundary (SEB) in the quasar 
mass-luminosity plane to compare different virial mass es- 
timators for quasar black hole masses. In particular, the 
Mgll estimator developed by McLure & Dunlop (2004) 
has been considered less reliable than the H/3-based es- 
timator of Vestergaard & Peterson (2006), with sev- 
eral proposals for corrections ■'Oiikeii £i Kollmcicr 20081 
|Risaliti, Young, fc Elvis 2009| IMarconi et al. 2009p . The 
quasar M — L plane indicates, surprisingly, that Mgll may 
be the most reliable virial mass estimator. A decline in the 
relative importance of non-virial motions at large radii may 
account for the differences in precision when using different 
emission lines. 

Surprisingly, using the adjacent continuum as a proxy 
for radius seems to be more precise than using time delays. 
If the statistical uncertainties of virial mass estimates are 
really smaller than previously believed, this would be a sub- 
stantial improvement. It means that we produce better black 
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hole mass determinations using one lower-resolution optical 
spectrum than from a time-series of high-resolution spectra. 

Moreover, many of the conclusions drawn from quasar 
mass functions, in Papers I and II, and in other work rely 
on the ability to segregate quasars into mass bins. With an 
uncertainty of 0.4 dcx, cross-contamination is a concern. An 
uncertainty of 0.15 dex allows us to divide a quasar sample 
into 2-3 times as many independent bins. Given the large 
systematic uncertainties in fitting bolometric luminosities 
to templates, it is possible that quasar mass functions will 
be more reliable than luminosity functions. 

The quasars that define the SEB are necessarily the 
closest to the Eddington limit of their cohort. It could be 
that this gives them a greater uniformity of properties, in- 
cluding more accurately virial motion in the broad-line re- 
gion, than quasars at lower Eddington rates. Lower Edding- 
ton rate quasars, such as those used in reverberation map- 
ping, might then have a wider dispersion in these same prop- 
erties, leading to the difference in observed spread. However, 
since the location SEB moves with redshift, something more 
complex than just the Eddington ratio would need to be re- 
sponsible for this possible greater uniformity in quasar prop- 
erties near the SEB. 

We also compared the SEB slopes using Mgll masses 
and H/3. Proposed corrections to Mgll masses comparing 
the two methods for the entire SDSS catalogue are quite se- 
vere, but these corrections might be due to uncertainties in 
the measurement of Mgll lines parameters. For the brightest 
objects at 0.4 < z < 0.8, the SEB produced by the Risal- 
iti, Young, & Elvis (2009) correction is a significantly worse 
match for the SEB produced by H/3 than using uncorrected 
Mgll-based virial masses. 

The authors would like to thank Bradley Peterson, Yue 
Shen, Michael Strauss, and Jonathan Trump for valuable 
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grant number G07-8136A. 
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